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Abstract 

In this work we rederive the Lamb-Retherford energy shift for an atomic electron in the 
presence of a thermal radiation. Using the Dalibard, Dupont-Roc and Cohen- Tannoudji 
(DDC) formalism, where physical observables are expressed as convolutions of suitable statis- 
tical functions, we construct the electromagnetic field propagator of Thermo Field Dynamics 
in the Coulomb gauge in order to investigate finite temperature effects on the atomic energy 
levels. In the same context, we also analyze the problem of the ground state stability. 
PACS numbers: ll.10.Wx, 32.80.-t. 



I. Introduction 



Since the 70's it has been argued [1] [2] that the physical interpretation of radiative phenomena, 
in particular the shift of atomic energy levels, rely upon different choices in the ordering of 
atomic and field operators in the interaction Hamiltonian. 

Almost two decades ago Dalibard, Dupont-Roc and Cohen- Tannoudji (DDC) [3] considered 
the interaction between a non-relativistic atomic electron and the quantized electromagnetic 
field, showing that the above mentioned arbitrariness can be removed by requiring the ob- 
servables variation rates to be Hermitian, if we want them to have a physical meaning. They 
generalized their procedure to the case of a small system S interacting with a large reservoir 
1Z (which may be in thermal equilibrium). This construct allowed them to separate the phys- 
ical processes in two categories, those where 1Z fluctuates and polarizes S (effects of reservoir 
fluctuations), and those where S polarizes 1Z (effects of self- reaction or radiation reaction). 

In the present work we are interested in analyzing the temperature effects in the context of 
DDC formalism, where the statistical functions, which are defined from two-point correlation 
functions, play a fundamental role. These functions enable us to obtain expressions, up to 
second order in perturbation theory, in terms of products of symmetrical correlation functions 
and susceptibilities [4] . The temperature implementation [3] can be made directly in such sta- 
tistical functions using the equipartition theorem, leading to a finite temperature description 
of the relevant phenomena. 

In an alternative way, we shall study the problem using Umezawa's formalism, known as 
Thermo Field Dynamics (TFD) [5]. In TFD, the quantum statistical average of a physical 
observable in a given ensemble is identified with its expectation value in a thermal vacuum. In 
this approach, temperature is introduced as an input in the eigenstates of the number operator 
associated to the quantized field. 

Our idea is to investigate the thermal propagator of the electromagnetic field in the 
Coulomb gauge in order to identify the symmetric correlation functions and susceptibilities 
of DDC formalism. This is the matter of section III, after a brief presentation of the main 
results of DDC construct in section II. In section IV we investigate the temperature depen- 
dence of the Lamb-Retherford energy shift of an atomic electron in the presence of a thermal 
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radiation field. In section V we calculate the variation rate of the mean atomic energy and 
discuss the stability of the ground state at finite temperature. In both cases, we are assuming 
that the whole system is in thermodynamic equilibrium. Finally, in section VI, we draw some 
conclusions. 

II. Radiation considered as a Reservoir 

In Dalibard, Dupont-Roc and Cohen- Tannoudji [3] construct, the interaction between an atom 
and the free electromagnetic field can be seen as the interaction of a microscopic system S 
with a large reservoir TZ, in the sense that TZ has many degrees of freedom and the correlation 
time among observables of TZ is small, allowing a perturbative treatment of the effect due to 
the coupling of S and TZ. 

Considering S an atom fixed at the origin of the coordinate system and TZ an homogeneous 
and isotropic broadband radiation field they addressed, among others, the problem of atomic 
radiative corrections as the Lamb shift and the dynamic AC Stark effect. In particular 1 , 
they showed that the shift in an atomic energy level (say a) caused by its interaction with the 
radiation field can be expressed as 

1 roo rj, j 

(SHsaf = -yr ^XkMW , (1) 

(5H S aT = ~ /_~ ^x'n^Cs^u) . (2) 

where Cs a (ou) (resp. Cr{uj)) and xsaiyi) (resp. Xr{oj)) are, respectively, the symmetric cor- 
relation and the even parity part of the susceptibility functions related to the atomic system 
(resp. reservoir) in frequency space. Their true meaning is well established in the context of 
DDC formalism which associates (1) to the reservoir fluctuation effects and (2) to the radiation 
reaction effects. 

Our main concern is the fact that such statistical functions are defined from two point 
functions of the dynamic operators involved in the interaction Hamiltonian Hj. In the case 

1 Here the total hamiltonian is given by H = Hs + Hr + Hi, where H$ describes the atomic system, Hr 
the radiation field and Hj is the interaction Hamiltonian. 
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we are interested in, the interaction Hamiltonian in the Coulomb gauge is given by 

H I = -(-)p.A(0), (3) 

where p is the momentum associated to the electron's motion and A(0) is the electromagnetic 
potential in the dipole approximation. The two point function for the spatial component of 
the field variable A^O, t) = A(t) (i = x, y, z) is given by 

^W^WO.^Wa + ^^VMA]-),, (4) 

where r = t' — t" and (,) A indicates an average on the reservoir state defined by a given 
statistical weight. As pointed before, the first term in (4) corresponds to the symmetric 
correlation function and the second is related to the linear susceptibility of the reservoir. The 
symmetric correlation function of the observable A(t), 

C R (T) = \({A(t>),A(t")} + ) A , (5) 

is real and tends to the ordinary correlation function in the classical limit. It gives a phys- 
ical description of the dynamics of fluctuations of the observable A(t). The other statistical 
function is the linear susceptibility Xr{ t )i which characterizes the reservoir response to an 
external perturbation, defined by 

Xr(t) = {o(r){[A(t%A(t")U) A 

= ^(r)Im^(-r), (6) 

where 6{r) is the step function. 

Since we are interested in analyzing the finite temperature dependency of (1) and (2), we 
postpone the calculation of the above statistical functions for the field components to the next 
section, where we shall employ the TFD formalism in order to obtain the finite temperature 
two point functions for the radiation field. 

Nevertheless, we restrict ourselves to the present action of the corresponding correlation 
and susceptibility functions in frequency space for the x component 2 of the atomic variable 
2 It can be shown that for i ^ j the statistical functions vanish [9]. 
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(ep/m), in the situation the atom is found in a given state \a) (with Hs\a) = E a \a)), namely 



e 2 



CfM = E —M Px \b)\\[5{u ab + + 5{u ab - u>)\, (7) 



h m 2 

-e 2 ,. ...... I _ I 



v + v- 

. U ab + LO LO ab - LO 



(8) 



= J2j^\(a\pM^(u; ab + u;)-6(u; ab -u;)}. (9) 

Expressions (8) and (9) are obtained by splitting the atomic susceptibility according to 
Xsaiyi) = x'sai 00 ) + iXsai 00 ) ' where each part characterizes, respectively, the response in phase 
and in quadrature at the frequency cu. In expression (8), V denotes the principal value. 

III. Thermal Correlation Functions and Susceptibilities 

In this section we study the thermal propagator of the electromagnetic field of Thermo Field 
Dynamics (TFD). Our idea is to define the statistical functions Cr(uj) and Xni 00 ) from the 
appropriated propagator of QED, implementing temperature at the beginning. We start by 
writing the quantized electromagnetic potential Ai(t) as 3 

A t (t) = A ( +\t) + A ( r\t), (10) 

where A^~\t) and A\ \t) are the components with positive and negative frequency, defined 
respectively as 

A ( + ) (t) = E« k e[(k)4e^, (11) 



1 - 

k,r 

with 



^ ( r ) W = E«k e [(k)«k t e^, (12) 



n \ 1/2 

In TFD we double the field degrees of freedom introducing the tilde conjugated of Ai(t) 
[5] [6]. Using the thermal doublet notation [6] [7], we obtain 



A ^)=(l;$ ) A l (t) = (A t (t), -A t (t)), (14) 



3 As in the last section, we assume that the atom is at rest at the origin of the coordinate system (r = 0) 
and that we are in the dipole approximation. 



where ( ) denotes the transposed and 



k,r 

A?\t)+A<r\t), (15) 



= E«ke[(k)(a;e^ + ^ e -^) 

k,r 

= A<+>(f)+A<r>( t ). (16) 



By construction, both fields Ai and A; are independent; the corresponding absorption and 
emission operators satisfy the algebra [6] 

[ «L ] = [ o k , 6 k f ] = 5 k ,k' <Jr )S - (17) 

At zero temperature, the vacuum state is given by the direct product |0) A <E> |0)^ = |0). 
Using (17), it follows that 

A { ? |0> = 0, i« |0> = 0. (18) 

In order to find the thermal propagator associated with the statistical functions, we must 
calculate the commutator 

[MVMnr = a^(t), (19) 

where fj,,v — 1,2 and i,j = x,y,z. The anti-diagonal components of the above quantity 
are identically zero when their expectation value in the |0) state is taken. The component 
/i — v — 1 can be written as 

Aj.l^AjWWlAjf-V), (20) 

where 

Af + \r) = [A^(t'),A^(t")}, (21) 
Af-\r) = [A^\t'),Af(t")]. (22) 
Now, using (11), (12), (17) and (18), we obtain, 

A^ (+) (r) = £ °£ <(k) ej(k) e-**\ (23) 



(24) 



k,r 



From (23) and (24), we can define two functional: 



A^ M) (r) = ^r)A^ (+) (r)+^(r)A^ ( - } (r) 



AH 

»j (ret) 



(+), 



r) + A^ 



(-), 



and 



A" ( i)(t) = A" 1 (r) - A 



ij (ret) 



'A'-'M. 



It can be easily shown that, at zero temperature, 



A^ (1) (r) =(0|{A(O,^-(O}|0>. 



(25) 



(26) 



(27) 



So, we should point out that the quantity g(r), defined before as a two time average of a given 
observable, is associated in the present case with the functional 



9a(r) = (0\M*)M? 



^. (1) (r) + ^.(r), 



(28) 



according to expression (4). 

By taking the Fourier transform of (25) and (26) we obtain, respectively, 



A^ (ret) M 



ii 



a*; (i)^) 



Va£ e[(k)e;(k) , , , 

^ ai <(k) e]{\i)Tt[8{u + u k ) + 5(u - 



(29) 
(30) 



Adopting the same procedure, we can extend the above calculation to the component /i — v — 
2. As a result, we have 



AgUH = E«keT(k)eJ(k) 



k,r 



ijj — — ie 



uj + — ie 



A?(i)H = -E«ke[(k)eJ(k>[% + Wk ) + ^-a; k )]. 



(31) 
(32) 



k,r 



We may write expressions (29) and (31) in a more compact notation, i.e., 

Aij (ret)M = V^e[(k)eKk) (- ; — 



(33) 
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and, in the same way, we write (30) and (32) as 



Aij = - ]T aJ&(k)eJ(k)7rT 8 [<S(w + co k ) + 5{u - u k )}, (34) 

k,r 

where, in the last two expressions, 

In TFD, it is known that the propagator at zero temperature is related to the one calculated 
in the thermal vacuum through a Bogoliubov transformation [8]. Applying this result to (33) 
and (34), we obtain, respectively, 

A&M = {B k l ^)A tl{ret) ^)B k ((3)r, (36) 



ij(ret) 
■ ij (1) 

where B k (f3) is give by 



AST&M = {B k \P)A lj{1) (uj)B k ((3)r\ (37) 



5 k (/3) = (l-n k ) 1/2 (^ _ f \- a f k J, (38) 
with a = 1/2, / k = exp[— fiuj k (3] and 

nk = = ^J37' ( 39 ) 

(P = 1/kT, where k is the Boltzmann constant and T the equilibrium temperature). 
The — u — 1 component of (36) is found to be 

+ i 7T [ 5(w k - w) - 5(cj k + cj) ] (l + 2n(wk)) }, (40) 

and, from (37), 

A " = E «k <(k) ej(k) tt [ % - u*) + 6{u + u k ) ] (1 + 2n{u k )). (41) 

k,r 

Now, relating (5) to (41) and (6) to (40), we are in position to define the thermal correlation 
function and susceptibilities, 

cgH = Al ia)H> ( 42 ) 
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and 



where, again, we split (43) as 



■ P, 



7? — _ — + p. 



(43) 

(44) 
(45) 



- £ o£ ej(k) e5(k) 

k,r 

Xl/M = "^ E ^<(k)e;(k)7r(l + 2n(a^)) [ <5(«k + ") - *(«k - ") ]■ (46) 

Choosing i = j = x and substituting the summation over modes by a polarization sum and 
an integral in k, we obtain 

j r^M 

C£M = ^— ^ / o dJ Hco'(n(co') + 1/2) [ 6(co> -u) + 6{cS + u) ] 

(47) 



37re c 3 

^%|KM) + i/2) 



37r£ c : 
1 



67re c 3 Vo 



doo'oo' 



V—^ — + V 



lo' — lo 



to' + CO 



x" x £{u) = 



U)M 



6%s c 3 Jo 

co (n(|w|) + l/2). 



dw'u' (2n(w') + 1) [ <J(w' + w) - - oo) 



3ne c 3 

IV. The Lamb-Retherford Shift via TFD 



(48) 



(49) 



We now apply the above results to the case of Lamb-Retherford shift and discuss the related 
thermal radiative effects. As mentioned in section I, this is done by substituting (8) and (47) 
into expression (1) which, according to [4], represent the desirable radiative correction, 



r W=^3EK«IpWI 2 x 



f^M 1 

x / duo'uo'iniuo') + -) 
Jo 2 



V 



1 



co ab + to' 



+ V- 



co ab - CO' 



(50) 



The atomic energy shift due to the field fluctuations appears as a sum of the effects (£# Sa ) fr ' of 
the "thermal photons", proportional to (n(u)), and that of the vacuum fluctuations (5Hs a Y v , 
corresponding to the u hu/2 by mode". This last term can be manipulated using the relations 



ru>M 



1 



/ co' duo' V- , 

JO CO ± UOq 



-r- i UJm i /nf u °\ 

lomT^oo hi h U{ 

uo K oo M y 



(51) 



/ J dJ 

Jo 



V—^ — +V- " 



= _2c ln^. (52) 



Hence, we obtain 



or 



« SH -r=<i^™m?^)^, (53) 

where a is the fine structure constant and hcK is the mean atomic excitation energy. Ex- 
pression (54) corresponds to the (pure) Lamb-Retherford shift as found in literature [10]. It 
is well known that its physical origin comes from the vacuum fluctuation of the radiation field 
(reservoir). The presence of k in (54) shows the quantum character of this effect, just as the 
vacuum fluctuation which gives rise to it. 

The contribution (5Hs a ) h ' proportional to correspond to a stimulated radiative 

correction due to the "thermal photons" . It resembles the AC Stark effect when the thermal 
radiation field is substituted by a quantized electromagnetic field. In the present context 
(5H S aY r ' corresponds to thermal radiative correction to the (pure) Lamb-Retherford shift and 
its effect vanishes as the temperature approach to zero. 

V. Energy Exchange 

In order to analyze the effects of the thermal reservoir on the stability of the atomic ground 
state, we now consider the energy exchange between a bound electron and the thermal radia- 
tion field using the results of section III. Following [4], the variation rate of the mean atomic 
energy when the system is in its ground state (say a) is given by 

j t (Hst = Y,{E b -E a )T a ^ (55) 

where T a ^b represents the transition rate between the ground state a and an excited state b 
due to the interaction with the reservoir. It is shown in reference [3] that (55) can be written 

as 

j t (Hsf a = Q^+Q-, (56) 
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where 

2/3 = /^C5( W )xS», (57) 
QJ3 = -/^^xSV)^). (58) 

The last two expressions have a clear meaning: (57) is associated with the energy absorption 
by the system when it is affected by reservoir fluctuations and (58) is related to the damping 
of the atomic motion caused by the reservoir. 

Using expressions (9) and (47) and taking into account the spatial components x, y and z 
of the electromagnetic potential, we find that (57) can be written as 

Qj = Q fr ' + Q fv 

= ^- E ^z[HM))+m, (59) 

b 

where 

_ e 2 |(a|p|6)|> afe | 

is the rate of spontaneous emission related to the transition between the levels b and a. 

The quantity Q" is calculated in the same way from expressions (7), (49) e (58). As a 
result, we find 

Q£ = Q fr ' + Q fv 

= -J2(E b -E a )rZ[{n(\u, ab \)) + l/2}. (61) 

b 

Substituting (59) and (61) in (56) we conclude that 

>>J = l (62) 

This result is what we must expect since the whole system is in thermal equilibrium at tem- 
perature T. In the present context, one can say that, in thermodynamic equilibrium, the 
bound electron reaches a new ground state which corresponds to the original one shifted by 
the amount (8H Sa )% For T = 0, the ground state stability still holds, since the effects of 
radiation reaction, Q£j, are cancelled by the effects of thermal vacuum fluctuations, Q%. 
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VI. Concluding Remarks 

In the present work we have used the structure of DDC construct to implement temperature 
effects via TFD. After a brief review of the main DDC results, we have investigated the 
propagators of the electromagnetic field in the context of TFD and derived the symmetric 
correlation functions and susceptibilities for the field variable A(0,t). Applying the results 
to the case of an atomic electron interacting with a thermal radiation field, we calculate the 
Lamb-Retherford energy shift and the corresponding corrections due to thermal photons. 

In the last section we have analyzed the energy exchange between the atomic electron and 
the thermal radiation field and concluded that, once the whole system is in thermodynamic 
equilibrium at a given temperature T, the stability of the ground state is maintained, even 
when T approaches to zero. 

We must point out that the original DDC formalism includes the case where the reservoir 
is a thermal radiation field. As remarked in [3], this is done by replacing the mean number of 
particles (n(|u;|)) by a Bose-Einsten distribution in the resulting statistical functions. How- 
ever, such procedure differs from ours in the sense that the detailed balance principle become 
meaningless in the context of TFD where the population dynamics between two given atomic 
states is not accessible. 

Finally, we mention that the applicability of TFD in the scope of DDC formalism is not 
restricted to the problem we have just revisited. Among the physical phenomena we intend to 
investigate in the near future are those related to the dissipative processes in quantum optics 
[11] [12]. 
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